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The Pauli operators (tensor products of Pauli matrices) 
provide a complete basis of operators on the Hilbert space of 
N qubits. We prove that the set of 4^ — 1 Pauli operators may 
be partitioned into 2^ + 1 distinct subsets, each consisting of 
2^ — 1 internally commuting observables. Furthermore, each 
such partitioning defines a unique choice of 2^ + 1 mutually 
unbiased basis sets in the iV-qubit Hilbert space. Examples 
for 2 and 3 qubit systems are discussed with emphasis on the 
nature and amount of entanglement that occurs within these 
basis sets. 

PACS numbers: 03.67.-a, 03.65. Ta, 03.65.Ud 



I. INTRODUCTION 

A pure quantum state of an A-qubit system is speci- 
fied by the eigenvalues of A independent commuting bi- 
nary observables ( "TV qubits carry A bits of information" 
In fact, a complete set of 2 N such states are so 
specified, each being associated with a binary number 
consisting of the A*" eigenvalues. Of the many alterna- 
tive choices of observable sets that define A-qubit basis 
sets, we are interested in those that are maximally in- 
compatible in the sense that a state producing precise 
measurement results in one set produces maximally ran- 
dom results in the other. 

The Pauli operators || provide an explicit realization 
of points raised above. First, they illustrate that al- 
though a greater number 2 N — 1 of observables simulta- 
neously take definite values, only A of these are required 
to define a pure state, and in fact these A generate all of 
the remaining compatible observables through multipli- 
cation. On the other hand, all A N — 1 Pauli operators are 
required in order to determine an arbitrary mixed state. 
In this connection, we shall show that the full set of these 
operators is exhausted in forming 2^ + 1 distinct subsets, 
each consisting of 2 N — 1 internally commuting observ- 
ables, and each defining its own unique eigenbasis. Both 
the observable sets and the corresponding basis sets are 
called mutually unbiased (J] (and in previous works the 
observable sets have also been called mutually comple- 
mentary Q)) because of the following physical property: 
If an A qubit system is prepared in a joint eigenstate of 
one such observable set, then it has a uniform probability 
distribution over the joint eigenstates of any of the other 
sets. It follows that all 2^(2^ — 1) observables outside the 
original (maximal) set of 2 N — 1 compatible observables 



will produce measurement results that are uniformly dis- 
tributed over all possibilities. Since the Pauli operators 
have binary spectra, it also follows that their dispersion 
is maximized. 

The equivalence of unbiasedness of basis sets and op- 
erator sets may be understood from the formal defini- 
tion as applied to basis sets, which may be summarized 
in general terms as follows: Let us denote basis sets by 
A = 1,2,..., and states within a basis by \A, a), with 
a = 1,2, ...,d (for the moment, we consider a Hilbert 
space with general dimension d, although our interest 
here is in d — 2 ). Two bases A and B are said to be 
mutually unbiased |^,^) if a system prepared in any el- 
ement of A (such as \A, a)) has a uniform probability 
distribution of being found in any element of B 



\(A,a\B,P)\< 



(1) 



where individual bases are understood to be orthonormal, 



{A,a\A,(3) = S a p. 



(2) 



Certainty of measurement outcomes for the operator set 
defining the \A, a)'s implies a uniform probability dis- 
tribution over states \B,[3), and this in turn implies a 
uniform probability distribution over all eigenvalue sets 
(and distinct measurement outcomes) of operators defin- 
ing the \B,/3)'s. 

A particular motivation for considering unbiased basis 
sets is that they provide for the most efficient determi- 
nation, using measurements alone, of a general (pure or 
impure) quantum state ||. In a d-dimensional Hilbert 
space, one needs d 2 — 1 real parameters to specify a gen- 
eral density matrix p, which must be hermitean and have 
Tr(p) = 1. Since measurements within a particular ba- 
sis set can yield only d — 1 independent probabilities, one 
needs d+ 1 distinct basis sets to provide the required total 
number of d 2 — 1 independent probabilities. Ivanovic || 
showed that the required number d + 1 of unbiased basis 
sets indeed exists if d is a prime number, and Wootters 
and Fields showed that it exists if d is any power of a 
prime number. Our proof is based upon this theorem of 
Wootters and Fields. 

The question of the existence and construction of unbi- 
ased basis sets is interesting not only from a fundamental 
point of view (e.g., in the formulation of "quantum me- 
chanics without probability amplitudes" 0, and in the 
information-theoretic formulation of quantum mechanics 
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0), but also as an important ingredient in quantum- 
information protocols (e.g., in the solution of "the mean 
king's problem" and in quantum cryptography 11). 
In particular, it was found recently that key distribu- 
tions based on higher-dimensional quantum systems with 
larger numbers of unbiased basis sets can have certain 
advantages over those based on qubits || . 

The present paper illustrates how the study of oper- 
ator relationships can provide a useful approach to the 
construction of unbiased basis sets of entangled as well 
as product character. The TV-qubit Hilbert space has 
dimension d — 2 N , and operators on this space (which 
include the density matrix) live in their own vector space 
of dimension 4^. The complete basis consisting of the 
Pauli operators || may be written as follows: Starting 
with the usual 2x2 Pauli matrices and the identity / 
that act on the spaces of individual qubits, 

o> = (<r x , <T y , a z , I), /i= (1, 2, 3, 4), (3) 

we write the 4^ tensor products (the Pauli operators and 
identity X) that act on the iV-qubit Hilbert space as 

JV 

°i = <T l(l,i) (T l{2,i)--- ,J fi(N,i) = LI 

k=l 

where k is the particle label and i distinguishes among 
the 4^ choices of the N subscripts fj,(k,i). This basis is 
orthonormal po| ; the inner product of two operators is 
defined as the trace of their product, 

Tr(OA) = II K<M)°2<W>) 
fc=i 

N 

= JJ 2 ^/*(fe,i)|*(fcj) = 2iV %, (5) 
fc=l 

where i = j means that fj,(k, i) = fj,(k,j) for every particle 
k. Like the individual Pauli matrices, each tensor prod- 
uct is self-inverse, Of — X, and apart from the identity 
(for which we reserve i — 4 N so that 4 n = J 1 / 2 ...I = 
X) they are all traceless and have eigenvalues ±1. 

The binary spectrum for each observable Oi permits its 
expression as a binary proposition: The two eigenvalues 
±1 of the observable Oi correspond to the values "true" 
or "false" of the proposition "The product of the spin 
projections cr^ (l i) a^ (2 A y-^ N ^ is +1." (If a particular 
az happens to be the identity, then no statement is made 
about the fcth qubit.) 

II. GENERAL RESULTS FOR N QUBITS 

We may now proceed to demonstrate the main formal 
points of the paper: First, that the set of 4^ — 1 Pauli 
operators (excluding the identity) may be partitioned 



into 2+1 subsets, each consisting of 2 — 1 internally 
commuting members, and second, that every such par- 
titioning defines a unique choice of unbiased basis sets 
(i.e., there is a one-to-one mapping from partitionings to 
choices of unbiased basis sets). 

The first part makes use of the proven existence of 
2 N + 1 unbiased basis sets ||. The projectors onto the 
unbiased basis states, 

P£ = \A,a){A,a\, (6) 

may be used to re-express the unbiasedness of bases A ^ 
B (Eq. [TJ) as 

Tr(P^) = 2- N . (7) 

and to define a set of operators O^ by their spectral 
decompositions 

Of = 5> Q P Q A (8) 

a=l 

We define e aa as a 2^ x 2^ matrix consisting of orthogo- 
nal row vectors, one of whose entries are all +l's, and the 
remaining of whose entries are equal numbers of +l's and 
— l's. There are exactly 2^ such orthogonal vectors, the 
components of each vector a being the eigenvalues of 0„ . 
One of these operators (say the a = 2 N th) is proportional 
to the identity, 0^ N = X. We include this to make Eq. 
(8) invertible, which will be useful later. The columns a 
label the joint eigenstates of the O^ (a = 1,2, 2 —1), 
and comprise the truth tables associated with the 2 N — 1 
corresponding propositions. This labeling is redundant; 
clearly an appropriate subset of just N rows may be used 
to construct iV-component column vectors that define all 
2 N joint eigenstates unambiguously as binary numbers. 
This reflects a property of the Pauli operators mentioned 
earlier. 

The above definition provides 2 N + 1 distinct sets (in- 
dexed by A), each set containing 2^ — 1 operators (fixed 
A and running index a = 1, 2^ — 1), after discarding 
the identity. Each of these operators has the spectrum 
±1 and is traceless, by construction. To show that they 
are unitarily equivalent to the Pauli operators, we need 
only demonstrate that they form an orthonormal set. For 
the case A ^ B, 

Tr(Of Of ) = £ e aa e bfi T!(P£P$) = 0, (9) 

where Eq. (Q) and the property ^ a e aa — were used; 
then for the case A = B, using Eqs. (||) and (||), 

Tr(Of Of) =J2z™s b pTr(P£p£) 

a, p 

= Y, E ac t £}> a = 2 N 8 ah . (10) 
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Finally, this orthonormal set of 4 — 1 traceless operators 
is completed by adding the identity, so indeed they have 
a representation in the form given byEq. (|). This shows 
that the Pauli operators may be partitioned accordingly. 

We now show the second part, namely, that any such 
partitioning of Pauli operators defines a unique choice 
of unbiased basis sets. Assuming such a partitioning, 
each subset (^4) of Pauli operators {0^,0^, ■■■,0^ N _ 1 } 
defines a unique basis of 2 N joint eigenstates \A, a), a = 
1, ...2^. Thus, each operator may be expanded as in 
Eq. (||), with e aa now defined as the eigenvalue of 
on the state \A, a), the lower index taking the values a = 
1,...,2 N — 1. The known spectrum of the 's dictates 
that each row of the s aa matrix must consist of an equal 
number of +l's and — l's, and the identity Tr(O^O^) = 
2 N S a b shows (note Eq. 10) that any two rows a and b 
are orthogonal. Thus, by appending an additional row 
(a = 2^) to the e aa matrix, we recover its previous form. 
The scaled matrix e aa / V2^ is orthogonal, and, therefore, 
we may invert Eq. (^) to yield the projection operators 

P^2- N Y.^0^2- N (l + Y.^). (11) 

a a 

In the second equality we write the identity contribution 
explicitly and delete the a = 2 N term from the sum, as 
denoted by the prime. 

We may now show that all of the basis sets are mutu- 
ally unbiased: Substituting Eq. ([ll]) into Eq. ([jj) yields 

Tr(P A P p B ) = 

2 -n + 4- N Y,Y, e aa SfMO^), (12) 

a b 

since terms linear in O a have vanishing trace. It follows 
immediately that if A and B refer to different basis sets, 
then Eq. ((?]) is satisfied. If A = B, then only the a = b 
term in the sum survives and 

Tr(P A P^ = 2- N (l + J2^ a s 0a ) 

a 

= 2~ N ^ £ aa£pa = $af3, (13) 
a 

where the orthogonality of £ aa j\f2~^ was used. This es- 
tablishes that the 2^ — 1 basis sets generated (uniquely) 
by the commuting subsets of Pauli operators are in fact 
unbiased. So there is a one-to-one correspondence be- 
tween partitionings of Pauli operators and choices of un- 
biased basis sets. 

III. EXAMPLES FOR TWO AND THREE QUBITS 

We now illustrate this correspondence for systems of 
two and three qubits. To develop notation, the operator 



subsets for the one qubit case consist of single elements, 
o~ x , a y , and a z . Corresponding basis sets are denoted by 
(x), (y), and (z), where each basis set consists of the two 
states "up" and "down" along the indicated axis. The in- 
dividual basis states are denoted by \n x ), \n y ), and \n z ), 
where n x = 1 or for spin "up" and "down" respectively. 
The inner products between any two states appearing in 
these basis sets obey Eqs. (1) and (2). Obviously, mea- 
surements by any of the above operators on an eigenstate 
of any other will produce perfectly random results (i.e., 
an average spin projection of zero). 

In the case of two qubits, the dimension of the Hilbert 
space is d = 4, so that five unbiased basis sets exist. 
Fig. 1 shows these together with the five corresponding 
operator sets, each consisting of three compatible oper- 
ators. Subscripts indicate three product bases, (zz)^, 
(xy)^, and (yx)^, whose individual states are denoted in 
the (zz),r case, for example, by \n 1 zl n 2 ). There are two 
Bell bases, (zx)b and (yz)si, and the states belonging to 
each of these may be written, respectively, as 

\nlnl;±) = -j=(\nl,n 2 x ) ±\nl,nl)), (14) 

W v ,nl\±i) = — (\n l v ,nl) ±i|n*,r^)), (15) 

where bars denote spin flips; ic, if n x = 1 or 0, then 
fix = or 1, respectively. Thus, the four individual basis 
states are explicitly enumerated [in Eq. (|i~5|), for exam- 
ple] by |l a ,l^;±i) and \l y ,O z ;±i) . The factor of i [as 
denoted by the subscript in the basis label (yz)Bi] is not 
arbitrary; its presence is dictated by the operators that 
define the basis, or equivalently by the requirement that 
the two Bell bases be mutually unbiased. It is a prop- 
erty of Bell bases that they can appear equally simple if 
other quantization axes are chosen in appropriate com- 
binations. For example, the Bell bases 4 and 5 in Fig. 1 
may be written as {yy)s and (xx)b, respectively. While 
these are simply different ways of writing the same basis 
sets, truely different alternatives involving all five basis 
sets exist for two qubits (see Refs. [[ll] and Jl^l). These 
alternatives also consist of three product bases and two 
Bell bases. 

The three-qubit Hilbert space has d = 8, and thus nine 
unbiased basis sets. One choice contains three product 
bases, (xyz)^, (yzx)^, and (zxy)^; and six bases con- 
sisting of maximally entangled states, {xxx)ci, {yyy)G, 
(zzz)a, (xzy)c, (yxz)c, and (zyx)c- The nine basis sets 
are listed in Fig. 2 and represented graphically in Fig. 3. 
The entangled basis sets are labeled by coordinate axes 
in which the states reduce to the familiar Greenberger- 
Horne-Zeilinger |l3| (GHZ) form. For example, all of the 
states belonging to the basis set (zyx)Q may be written 
as 

\nl,n 2 y ,n 3 x ;±) = -^=(\nl,n 2 y ,nl) ± \n l z ,n 2 y ,nl))] (16) 
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FIG. 1. Five unbiased bases sets and corresponding Pauli 
operator sets. Each operator set consists of three commuting 
members, any two of which determine the corresponding basis 
set as their joint eigenbasis. 
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FIG. 2. Listing of nine unbiased basis sets and correspond- 
ing operator sets, each consisting of seven commuting mem- 
bers. Particular subsets of three determine the corresponding 
basis sets completely. 



these would require more complicated expressions if re- 
ferred to other coordinate axes. Fig. 2 lists the seven- 
member operator sets that correspond uniquely to each 
basis set. As in the two-qubit case, the operators involv- 
ing only a single Pauli matrix are exhausted within the 
product basis sets. 

It is striking that in the progression from one to two to 
three qubits, the number of totally entangled bases can 
grow from none to two to six, while the number of prod- 
uct bases remains fixed at three. It is easy to convince 
onesself that the maximum number of product bases re- 
mains fixed at three for all numbers N of qubits. 

To show that the structure is more flexible with three 
qubits than with two, we describe now a different choice 
of unbiased basis sets for three qubits, one that cannot 
be obtained from the previous choice by local unitary 
transformations. In this choice (Fig. 4), there are no 




3 product bases 6 GHZ-type bases 

FIG. 3. Schematic of unbiased basis sets listed in Fig. 2 - 
three product and six GHZ bases. The three particles in the 
circle are maximally entangled (in a GHZ state). 
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FIG. 4. Same structure as in Fig. 3, but all bases are 
partially entangled. In each of three groups a different particle 
must be singled out as unentangled. 
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FIG. 5. Schematic of unbiased basis sets listed in Fig. 4, 
all partially entangled. One particle is unentangled with the 
other two, which are totally entangled in Bell states. 



product states, and no states with three-particle entan- 
glement. Every basis consists of states that are prod- 
ucts of one-particle states with Bell states; one particle 
is unentangled while the other two are totally entangled, 
as depicted in Fig. 5. The basis sets form groupings of 
three: (x 1 )(yz) B i, (y 1 )(zx) B , (z 1 )(xy) B , then (x 2 )(xy) B , 
(y 2 ){v z )Bi, (z 2 )(zx) B , and finally (x 3 )(zx) B , (y 3 )(xy) B , 
(z 3 )(yz) B i, in which a different particle is factored out 
within each group. Coordinate axes are permuted within 
each group, but not from group to group. Factors of 
% appear once within each group. 

Note that within each grouping, we find three unbi- 
ased Bell-type bases - a feature that was not seen in the 
two-qubit system. Indeed, if one were to begin with three 
unbiased Bell bases in a two-qubit system, one could then 
not find two additional basis sets. This can be understood 
in terms of the operator decomposition: The nine opera- 
tors exhausted by three Bell bases do not leave six opera- 
tors that are decomposable into two commuting subsets. 

We also note that the choices of unbiased basis sets 
given for the three-qubit case are not obtained from the 
algorithmic construction given in Ref. [§. Wootters has 
pointed out jl4|] that this construction produces another 
choice which consists of two product, four GHZ, and 
three "product-Bell" bases. In the case of two qubits, 
not surprisingly, the same construction produces three 
product and two Bell bases [H. 
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IV. CONCLUSIONS 

The two- and three-qubit cases illustrate general points 
made at the beginning. First, with regard to state prepa- 
ration, one can see that N observables suffice to define 
any of the listed basis states completely, representing 
these states as binary numbers. In the two-qubit case, 
any two of the three compatible observables within a 
subset may be chosen. In the three-qubit case, there 
are many choices of three observables that suffice - for 
example, the first three listed within each subset. In the 
iV-qubit case, we introduced an e aa matrix in which an 
appropriate choice of N rows (representing N operators) 
describe all 2 N basis states as binary nunbers. 

Second, with regard to the determination of a general, 
possibly mixed state, recall that 4^ — 1 real parameters 
(15 for two qubits and 63 for three qubits) are required 
to specify the A-qubit density matrix completely j|,|]]. 
And exactly this number is provided, either by the ex- 
pectation values of the operators themselves, or by all 
the independent probabilities associated with the unbi- 
ased basis states. As we have proven for the general case, 
the 4^ — 1 Pauli operators can be partitioned into 2^ + 1 
subsets, each consisting of 2 N — 1 internally commuting 
observables. The set of all such partitionings has a one- 
to-one correspondence with choices of 2 N + 1 unbiased 
basis sets in the A-qubit Hilbert space. There are many 
such choices, and for N > 2 the entanglement may be 
distributed over basis sets in many different ways. The 
maximum number of product bases is fixed at three for 
any N. 

The correspondence between basis sets and observables 
makes it possible to regard all Pauli operator subsets 
within a given partitioning as being mutually unbiased: 
If the system is prepared in a joint eigenstate of one ob- 
servable set, then it has a uniform probability distribu- 
tion over the joint eigenstates of any other observable set 
in the partitioning. As a result, all observables outside 
the original maximal commuting subset yield minimal 
information - measurement outcomes are uniformly dis- 
tributed over all possibilities. 

The concept of unbiasedness between observable sets 
extends the idea of complementarity of two individual 
observables that fail to commute. Clearly two such ob- 
servables must always belong to different mutually un- 
biased subsets within any partitioning. However, as the 
two- and three-qubit examples show, two commuting ob- 
servables may belong to the same or to different unbiased 
subsets. Their compatibility is dependent upon the par- 
titioning. 

Note added. After this work was completed, an e-print 
Jl2j appeared reporting work which is related to this 
work, but complementary in several respects. Ref. JO] 
obtained a general relationship between complete bases 
of unitary operators (beloning to the general Pauli group) 



and unbiased basis sets, for any power-of-prime dimen- 
sion. In this paper, we considered the many-qubit case. 
We expanded upon the physical interpretation of the con- 
cept of complementarity. We showed that many alter- 
native paritionings are possible and, most importantly, 
entanglement is distributed among unbiased basis sets in 
a partition-dependent manner for N > 2. 
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